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NOISE GENERATOR 


DECUS Prog ram Library Write-up 


DECUS No. 8-300 


ABSTRACT 

The program, which is written as a subroutine, creates 
a pseudo-random voltage with a gaussian probability den¬ 
sity function. Also the appropriate binary noise is available. 
The bandwidth can be selected by the programmer. 

REQUIREMENTS 

Storage 

The program requires one page of the memory. 

Equipment 

Basic PDP8 with 12 -bit Digital Output (Device Code No. 30) 
and one Digital-to-Anlog Converter (Device Code No. 55). 

USAGE 


Loading 

The program is written as a subroutine and ends with a 
PAUSE Statement. It must be assembled with the main program 

Calling Sequence 

Initially call the program with "JMS NOIINI". This part sets 
the timer and the register REGA is loaded with the content 
of the switoh register. Return is to the next instruction. 
Normally call the generator with "JMS NOIGEN". The number 
of callings per second determines the bandwidth of the 
noise voltage. 

Switch Setting 

The oontent of the switoh register is strobed into the 
accumulator, if the initial subroutine is called. It de¬ 
termines the presetting of register A. In the subroutine 
NOIGEN the switch register is never more used. 

Start Up 
None. 










T a 8.388.607 « 4- . 

J T 

If the subroutine is called once in 20 ms, the bit pattern 
is repeated not before 46 hours. 

The power spectrum of this binary noise is equal to the 
sin 2 x/x 2 function with the first zero at f T . 

The next part of the program uses the binary noise to create 
analog noise with a ganssian amplitude distribution. For this 
the binary noise is filtered by an ideal rectangular low 
pass filter with a bandwidth of about i/15 of f T . The pulse 
response of such an ideal filter is 

r ^ Sin 2r(t-t„)/f g 

(t) a T(t-t 0 ;/fg 

The response of any inputsignal to this filter can be 
obtained by using 

GO 

y (t) -/•(«) * * (t-u) du 
or by replacing the integral by a discret adding 



The part x(t-nAt) means that the input voltage must be ^ 

delayed by predetermined values. Since in the program x ( 
is a binary pattern and comes from a shift register with* 

36 bits the above equation can be performed on a simple 
way. The coefficients a n determine the values of (g, v. 

The limited length of the shift register do not allow to 
integrate to the infinite but only to n s 36. So the 
function sinx/x is only used between x « -3ir and x =* +3f 
The remaining part is seperated into 36 intervalls and 
the appropriate values are stored in a table. The values 
are calculated by an appropriate FOCAL Program which uses 
the equation: 
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Miscellaneous 


The following labels are used in this subroutine 


ANALOG 

NOIDUR 

BITCTR 

NOIGEN 

BIT1IIG 

NOIINI 

EBIT23 

RAND 

EBIT5 

RAOFFS 

EXOR 

REGA 

EXRES 

REGB 

LOADRA 

REGC 

Mi 2 

TABBEG 

NOICTR 

TABENT 


RESTRICTIONS 

None. 

DESCRIPTION 

Discussion 

The program makes use of a 36 bit shift register, con¬ 
taining three words named REGA, REGB and REGC. Each time 
the core, named NOICTR, overflows the full register is 
shifted one bit to the right. The input of the register 
is determined by the subroutine EXOR, which performs a 
logical exclusive-or-gate from bit 5 and 23 of the shift- 
register. One part of the siiiftregister is transferred 
to the Digital Output Register No. 30 . and gives a bi¬ 
nary noise with the two levels zero and one. This voltage 
has a cycletime which depends on the callings per seconds. 
The period of a whole cycle is calculated by 

T a (2 n - 1) * 4- 

i T 

with n number of bits of the shiftregister 

frp calling frequency of the noise generator. 

Since the exclusive—or—operation is performed between bit 5 
and bit 23 the number n becomes 23 and 






with 


n c I*2j •*«f 37 

The filter caw easily be changed by loading the table 
with other coefficients. The program masks the Binary 
pattern of the shiftregister with the table and adds 
the appropriate values to form the output voltage. To 
get a mean of zero volt an offeet voltage is added before 
converting the voltage by the 10 bit D/A converter. 

The amplitude of the analog noise is 1,15 V RMg and the 
maximal possible value is 4,87 V. So the crest factor 
becomes 4,2. 

Example 

The example in the appendix shows the binary noise and the 
analog noise. Also the correlation function and the power 
spectrum of both are calculated. Last not least the 
probability density function of the analog noise is shown. 

PROGRAM 


Program listing 

A listing of the source is given in the appendix 






EXAMPLE 
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Output voltages of the Noise Generator 



























































EXAMPLE 



6 


ainplefrequency 










EXAMPLE 
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EXAMPLE 
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Sample frequency 25 cps, N = 3000) Class \vidth 














0001 


/ DFVLf -Hivh.-L-LA 

/ USE fi9 - £•)3 
/ NOISE GENERATO i 
/ 81 JULY 19 69 
/ PDP8, BASIC CON FI PUPATION 

/ 

/ CALL SDUOIHINE BY JM S MO I GEN. 

/ ' FT URN IS TO FOLLOWING^ IN ST RDCTI0 N . 

/ IN t "i I ALL { CALL JM S NO IINI. 

/ 01 ! Ti'« *T IS VIA REGISTER #30 FOE. BINARY NOISE 
/ AN R W| (\ |}/ft- CONVERTER #55 FOR ANALOG NOISE. 

/ 


) I IN I , O' 

1.1. NO I Di IF: 
DGA NO I Cl « 

L A S 
STL 
hA R 

DC A Ah (ft 
CCA 1 FGh 
CCA REGC 
ON; I NOUNI 

/ 

NQIC-I.N, i 

I 7 NO I Cl 
JMR I MOIGF.N 
JMS f'.xo, 

lAi) REG A 

; at. 

L>CA EE GA 

1 1 \J r t t-r> 
h A h 

CCA HE GH 
1 ftb KEG i 
CAN 
630? 

DC A f. F GL 
JMS ANALOG 
TAD NO I DUE 
DC A NO I GTE: 
JMr I NO I PEN 

0 
0 
P 
0 

- 1 

0 


CCA EXCEL 
TAD REGA 
AND EBIT5 
SZA CLA 
ISZ EX HE S 
TAD EE GB 


/ 

EE FA, 
EEGB, 

K E GO, 
NOICTE, 
NO I DUE, 
/ 

F AO i , 


/ I NI 11 AL I Z E 

/ SET COUNTER 
/ GET CONTENT OF SR 
/ SET LX b 

/ SHIFT ONE BIT TO THE RIGHT 
/ S 10 KE IN RE GI STER A 
/ CLEAR REGISTER R AND C 


/ NOISE GENERATOR 
/ TINE TO CREATE VALUE ? 

/ NO 

/ EXL. Of OPERATION OF BIT S AND 23 
/ RESULT IS IN LKB 

/ SHIFT REGISTER a 


/ SHIFT REGISTER U 


/ SHIFT REGISTER C 

/ L0 AP i E GI Hi ER ii 30 Ivl TH LIN. NOISE 

/ GET ANALOG NOISE 

/ RESET COUNTER 
/ EX I 1 


/ REGISTERS 
/ TIMING COUNTER 
/ COUNTER BEGIN 


/ EXCLUSIVE OR OPERATION 
/OF BIT 5 AND BIT 23 
/ RESULT IN LINKBIT 
/ CLEAR 

/ GET REGISTER A 
/ MASK BIT 4 
/ BI T HI GH ? 

/ YES, LOAD EXRES 
/ GET REGISTER B 
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0002 


/ 

EXRES, 
EBI T5, 
EBIT23, 
/ 

ANALOG, 


/ 

LOADRA, 


H I THI G, 


/ 

M 1 2, 
RAND, 

BI TOTE, 
I A BENT, 
/ 

DECIMAL 
h AOFFS, 
1 ARRI' G, 


AND 

FBIT23 

/ 

MASK BIT 10 

SZA 

CLA 

/ 

BI T HIGH ? 

I SZ 

EX RES 

/ 

YES, EXRES+1 

TAD 

EX RES 

/ 

RESULT OF OPERATION IN BIT 11 

RAH 


/ 

SHIFT IN LKB 

CL A 




JMP 

I LXOE 

/ 

EX I T 

0 


/ 

BIT 11 CONTAINS RESULT 

200 


/ 

MASK FOR BIT 5 

2 


/ 

MASK FOR BIT 23 

0 


/ 

SET ANALOG NOISE VOLTAGE 



/ 

BY DIGITAL FILTERING 

TAD 

TAB ENT 



DC A 

BITCTR 

/ 

SET TABLE ENTRY 

DC A 

RAND 

/ 

CLEAR OUTPUT 

TAD 

M 12 



DCA 

NOICTR 

/ 

SET COUNTER TO FULL WORD 

TAD 

REGA 

/ 

GET REGISTER A 

JMS 

LOADFA 

/ 

LOAD OUTPUT 

TAD 

M 12 



DCA 

NOICTR 

/ 

RESET COUNTER 

TAD 

REGB 

/ 

GET REGISTER B 

JMS 

LOADRA 

/ 

LOAD OUTPUT 

TAD 

M 12 



DCA 

NOICTR 

/ 

RESET COUNTER 

TAD 

RE GO 

/ 

GET REGISTER C 

JMS 

LOADRA 

/ 

LOAD OUTPUT 

TAD 

RAND 

/ 

GET OUTPUT 

TAD 

RAOFFS 

/ 

ADD OFFSET VOLTAGE 

6552 

/ 

CONVERT TO ANALOG BY IV A #55 

CLA 




JMP 

I ANALOG 

/ 

EXIT 

0 


/ 

LOAD OUTPUT WITH APPROPRIATE 



/ 

VALUE FROM TABLE C SIN X / X ) 

CLL 


/ 

CLEAR LINKBIT 

RAL 


/ 

SHIFT BITS IN LKB 

SZL 


/ 

CHECK IF BIT IS ZERO 

JMP 

BITHIG 

/ 

FOUND BIT HIGH 

I SZ 

BITCTR 

/ 

IN CREM EN T TABL EPO INTER 

I SZ 

NOICTR 

/ 

HAVE ALL 12 BITS ? 

JMP 

LOADRA+1 

/ 

NO 

JMP 

I LOADRA 

/ 

YES, EXIT 

DCA 

EX RES 

/ 

FOUND BIT HIGH, SAVE AC 

TAD 

I BITCTR 

/ 

GET VALUE FROM TABLE 

TAD 

RAND 

/ 

ADD TO OUTPUT APPROPRIATE VALUE 

DCA 

RAND 

/ 

STORE- 

TAD 

EXRES 

/ 

GET saved ac 

JMP 

LOADRA+5 



- 1 4 


/ 

BITS PER WORD 

0 


/ 

STORAGE FOR OUTPUT 

0 


/ 

POINTER TO TABLE 

TABBED 

/ 

TABLE ENTRY 

-319 

/ 

OFFSET 

0 


/ 

TABLE FOR CFSIN X)/X 

6 


/ 

FOR X=PI/6.*(N- 19) 


in 










0003 


10 
1 3 
12 

7 

0 

-9 

- 1 7 
-21 
-21 

- 14 

0 

19 

41 

64 

8 3 
96 
100 
96 
8 3 
64 
41 
19 

0 

- 1 4 
-21 
-21 

- 1 7 
-9 
0 

7 

12 
1 3 
10 

6 

/ 

OCTAL 

FA USfv 





